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Abstract: 

The Hawking flux from a black hole, (at least as seen from asymptotic infinity), is 
extremely sparse and thin, with the average time between emission of the successive 
Hawking quanta being many times larger than the natural timescale set by the energies 
of the emitted quanta. While this result has been known for over 30 years, it has 
largely been forgotten, possibly because many subsequent authors focussed mainly on 
the late-time high-temperature regime. We shall instead focus on the early-stage low- 
temperature regime, and shall both quantify and significantly extend these observations 
in a number of different ways. In particular we shall confront numerical estimates with 
semi-analytic approximations based on a naive Planck spectrum. 

First we shall identify several natural dimensionless figures of merit, and thereby 
compare the mean time between emission of successive Hawking quanta to several 
distinct but quite natural timescales that can be associated with the emitted quanta, 
demonstrating that very large ratios are typical for emission of massless quanta from 
a Schwarzschild black hole. Furthermore these ratios are independent of the mass of 
the black hole as it slowly evolves. We shall then show that the situation for the more 
general Reissner-Nordstrom and generic “dirty” black holes is even worse, at least as 
long as the surrounding matter satisfies some suitable energy conditions. The situation 
for the Kerr and Kerr-Newman black holes (or even for charged particle emission from 
a Reissner-Nordstrom black hole) is considerably trickier, and depends on a careful 
accounting of the super-radiant modes. 

Overall, the Hawking quanta are seen to be dribbling out of the black hole one at a 
time, in an extremely slow cascade of 2-body decays. Among other things, this implies 
that the Hawking flux is subject to “shot noise”. Observationally, the Planck spectrum 
of the Hawking flux can only be determined by collecting and integrating data over a 
very long timescale. We conclude by connecting these points back to various kinematic 
aspects of the Hawking evaporation process. 
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1 Introduction 


It is (or should be) well-known that the asymptotic Hawking flux from a black hole is 
extremely sparse and extremely thin. The interstitial gap, the average time between 
emission of successive Hawking quanta, is many times larger than the natural timescale 
set by the energies of the emitted quanta themselves. This result was established over 
30 years ago [1-5], but has largely been forgotten, quite possibly because many of the 
subsequent authors focussed mainly on the late-time high-temperature regime in the 
final stages of the evaporation process [6-11]. The early-stage low-temperature regime 
has recently been reconsidered by van Putten [12, 13]. 

We also shall focus on this early-stage low-temperature regime, and shall develop 
several simple semi-analytic estimates based on assuming an exact Planck spectrum, 
which is not the full story but is sufficient to give tolerable estimates, at least for 
spin-zero bosons. We shall compare these semi-analytic estimates with the numerical 
estimates of Page [1-4]. (These numerical estimates in particular include the effects of 
greybody factors.) Overall, and carefully separating out super-radiant contributions, 
(which in Hod’s article [14] were simply lumped in with the Hawking effect), the sparsity 
of the Hawking flux is seen to persist throughout the entire evaporation process. 

We shall carefully define and justify several natural dimensionless figures of merit, 
which are suitable for comparing the mean time between emission of successive Hawking 
quanta with several natural timescales that can be associated with the emitted quanta. 
We shall first focus on non-super-radiant situations: We shall see that large ratios 
(very much greater than unity) are typical for emission of massless quanta from a 
Schwarzschild black hole. Furthermore these ratios are independent of the mass of 
the Schwarzschild black hole as it slowly evolves; certainly for as long as the Hawking 
temperature is well below the QCD scale. We shall then show that the situation for 
the more general Reissner-Nordstrom and generic “dirty” black holes is even worse, at 
least as long as the surrounding matter satisfies some suitable energy conditions, and 
as long as one is looking at the emission of uncharged Hawking quanta. We shall also 
consider the effects of particle rest mass on the emitted Hawking quanta. 

The onset of super-radiance considerably complicates the discussion. Super-radiance 
can (sometimes) occur for charged particle emission from Reissner-Nordstrom black 
holes, and for a range of emitted quanta from rotating black holes. In particular the 
situation for the Kerr and Kerr-Newman black holes is quite tricky, and depends on a 
careful accounting of the super-radiant modes. Certainly the quanta emitted in super¬ 
radiant modes are “quantum vacuum radiation”, but whether or not one chooses to 
call them Hawking quanta is more problematic. 
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Overall, throughout the entire history of the Hawking evaporation process, the 
(non-super-radiant) Hawking quanta will be seen to be dribbling out of the black hole 
one at a time, in an extremely slow cascade. Among other things, this implies the 
presence of “shot noise” in the Hawking flux. Observationally, the Planck-shaped 
spectrum of the Hawking flux can only be determined by collecting and integrating 
data over a very long time. We shall conclude by connecting these points back to 
various kinematic aspects of the Hawking evaporation process, which is now seen to 
resemble a cascading chain of 2-body decay processes. 

2 Strategy 

We shall compare and contrast two approaches: 

• As a zeroth-order approximation it should be perfectly adequate to treat the 
Hawking flux as though it is simply Planck spectrum blackbody emission at the 
Hawking temperature. While we know that a more careful treatment should at 
the very least include greybody, phase-space, and adiabaticity effects [15], (see 
also references [16, 17] and [18, 19]), nevertheless a zeroth order approximation 
using a blackbody spectrum should be quite sufficient to set the scale (if not the 
precise details) for the relevant issues we wish to consider. 

• At the next order of approximation, the most significant limitation on naively 
treating the Hawking flux as Planckian blackbody emission emission comes from 
the greybody factors. Page [1-4] resolves the Hawking flux into spin-dependent 
angular-momentum modes, and calculates various quantities of the form 

(Q) = Y. f Tsim{u) (n)u Q(oj) dw. (2.1) 

im 

Here (n) u is a completely standard bosonic/fermionic occupation number, while 
the Tstmiu) are spin-dependent greybody factors, estimated by numerically solv¬ 
ing the appropriate Regge-Wheeler/Zerilli equation for the radial waveform, this 
all being followed by a numerical integration over frequencies. 

Including the complications due to the integral over greybody factors does change 
the numerical value of our estimates, sometimes quite drastically, but does not 
change the qualitative nature of our results. (The much smaller effects due to 
adiabaticity and phase space constraints will remain negligible up to the final 
stages of the Hawking cascade when the black hole mass has shrunk down to the 
Planck regime.) 
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Between them, these two approaches give a good qualitative and quantitative handle on 
the sparsity of the Hawking flux. The blackbody emission approximation will be seen 
to work best for spin-zero, with higher spins seeing extra suppression (and increased 
sparsity) due to the angular momentum barrier. 


3 Flat space preliminaries 


The differential number flux, (quanta) / (time ), (of massless bosonic quanta emitted by 
a black body of temperature T, infinitesimal surface area d A, and surface normal h), 
into a wave-number range d 3 k is (in flat space) given by the utterly standard statistical 
mechanics result: 


dT = 


9 


c (k ■ n) 


d 3 k dA 


(3.1) 


(27r) 3 exp (hck/ksT) — 1 
Here g is the spin degeneracy factor; which is 1 for scalar bosons and 2 for massless 
bosons with non-zero spin. Then integrating over azimuthal directions, 


dT 


g c 27 rk 2 cos 6 

(27t) 3 exp (hck/k B T) — 1 


smOdOdk dA. 


(3.2) 


Integrating over the remaining angle, 6 e (0,7r/2), we see 1 


d r 


g ck' 2 

87T 2 exp (hck/k B T) — 1 


dk dA. 


(3.3) 


The wave-number integral can easily be performed, so that for an object of finite surface 
area A the total emitted number flux (see e.g. Schwabl [21]) is: 


g_m kyr 3 

47T 2 h 3 C 2 


(3.4) 


The reciprocal of this quantity, r gap = 1/T, is the average time interval between the 
emission of successive quanta. 

In counterpoint, the peak in the number spectrum occurs where k 2 / ( e hck / k B T _ l) 
is maximized, that is, at 


^ peak number 


ck 


peak number 


k B T 

h 


(2 + W(-2e~ 2 )) . 


(3.5) 


Here W(x) is the Lambert IH-function, defined by W(x)e w ^ = x. (See, for instance, 
references [22-24]; the presence of the Lambert W function in this calculation is not 

1 This angular integral is ultimately responsible for the | in the relationship between Stefan’s 

constant and the radiation constant: cr = \ac. See e.g Roberts [20]. 
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a “deep” result; it appears for exactly the same reason that the Lambert bb-function 
appears in the constant characterizing Wien’s displacement law.) The quanta emitted 
at this peak can only be temporally localized to within a few oscillation periods, so it 

is Safe to take ^localization 1 /7peaknumber number clS a good estimate of the 

time required for each individual quantum to be emitted. 2 Let us now define the 
dimensionless figure of merit 


^peaknumber _ 7T (2 + W ( 2e ' )) H C 


V peak number 


gap 


^localization 


9(13) 


k 2 B T 2 A' 


(3.6) 


In terms of the so-called “thermal wavelength”, A t h e rmai = 27T hc/(k B T), this is 


V peak number 


(2 + W(-2e- 2 )) Al h , rm „ 
4ir g C(3) A 


(3.7) 


If instead we consider the peak in the energy flux, rather than the peak in the number 
flux, then the only change is that now the factor (2 + bb(—2e -2 )) —* (3 + bb(—3e -3 )), 
and 


^ peak energy 


ck 


peak energy 


k B T 

h 


(3 +W(-3e- 3 )) . 


(3.8) 


Then we have 


(3 + W(~ 3e" 3 )) A thermal 

t? peak energy ^ ^ 


(3.9) 


Similarly, we could use the average frequency to set the localisation timescale 


/ x _ Jck(dT/dk)dk _ tt 4 k B T 
{ ' /(dT/dfc)dfc 30C(3) h ' 1 ' ’ 

The net result is that (2 + bb(—2e -2 )) —* 7t 4 /(30 C(3)), and so 

2 \ 2 

n —_ - _ A thermal (o i i A 

'/average energy — g £(3) 2 A ' V J A 


Now consider something more subtle; let us divide the spectrum into “wave-number 


bins” and take 


V binned 


1 

7 2t r dT 

J ck d k 



(3.12) 


2 Originally we had used Ti oca i iza tion = l/w pea k number ■ In counterpoint, Hod, see reference [14], 

prefers to use ri ocalizat ion = l/v peak number = 2-7T/oipeak number • This numerical factor does not qualita¬ 
tively change our results, but in the interests of being as conservative as possible we shall include the 

2tt. 
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This quantity effectively calculates the decay rate into wave-number bins, of width d k, 
centred on u> = ck , compares this with the frequency v = cu/(27r) = ck/(2 tt), and then 
sums over all bins. A brief calculation yields 


24 h 2 c 2 

Vanned ~ ^ 


9 A \ 2 
^ A thermal 

8n 3 g A 


(3.13) 


All of these sparsity estimates (in flat Minkowski space for now) take the form 

\2 

= (dimensionless number) the ™ al . (3-14) 

Let us now introduce key aspects of black hole physics, adapting the discussion above 
to see how far we can get. 


4 Non-super-radiant situations 

4.1 Schwarzschild black holes 

Under normal laboratory (and astronomical) conditions one is dealing with emitters 
whose surface area is extremely large in terms of the thermal wavelength, so in those 
situations r] • Cl. However, this is exactly what fails for a Schwarzschild black hole. 

• First T —>■ Th, and for the Hawking temperature we have 

fxc 

ksTn = ^ thermal = 87T 2 r H . (4.1) 

47T rn 

Note that the thermal wavelength is a factor 87T 2 ps 78.95 ~ 80 times larger than 
the Schwarzschild radius. 

• Second, a subtlety arises here as to which “area” to use. Naively one might use 

A —>■ Ah = 47T rjj, but this corresponds to a cross section of \Ah = nr 2 H , which 
is really only appropriate for some particle species in the low-frequency limit. At 
high frequencies (the ray optics limit) the cross section is universally given by 
T 7rr rr = 25]. This is enhanced by a factor of and implies that, (to 

smoothly match high frequency results), we should set A —)■ A e ff ective = ^ A H = 
27nr 2 H . 3 

3 This numerical factor does not qualitatively change our results, but in the interests of being as 
conservative as possible we shall include the . 
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4.1.1 Massless Bosons 

With these substitutions, for a Schwarzschild black hole we have 


A 

A 


2 

thermal 

effective 


647T 3 


« 73.49635955... > 1, 


(4.2) 


which is certainly much larger than unity. Consequently 

_ 32vr 2 (2 + W(-2e~ 2 )) 15.50768123... 

V peak number 

As promised, the gap between successive Hawking quanta is on average much larger 
than the natural timescale associated with each individual emitted quantum. Note that 
this is a physical situation where the relevant dimensionless constant is not “natural” - 
physically it is not of order unity, though mathematically one would quite legitimately 
still say the factor is 0(1)] neither is it zero nor infinity — the various numerical factors 
are important in determining the order of magnitude. Physically even more important 
is the fact that the mass of the Schwarzschild black hole drops out of the calculation, 
so that this calculation will be relevant as long as the dominant Hawking emission is 
into massless quanta. 4 This certainly holds throughout early stages of the evaporation 
where ksTn -C m e c 2 , and will plausibly remain relevant until IzbTh < Aqcd- 

Similar calculations apply for the other options we had considered for the localiza¬ 
tion timescale. Still working with the Schwarzschild black hole, we see that: 


27gC(3) 


i. 




• If we consider the peak in the energy flux, rather than the peak in the number 
flux, then 


V peak energy 


32tt 2 (3 + W(—3e“ 3 )) _ 27.45564528... 
27gC(3) ~g 


» 1 . 


(4.4) 


• If we consider the average frequency then 

_ 26.28537289... 

V average energy 


(4.5) 


• For the binned version of the p parameter we have 

14.22222222... , 

binned — • (4.6) 


4 We shall soon see that for other more general black holes the Schwarzschild result provides a 
mass-independent lower bound. 
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However one chooses the precise details to set up the calculation, and whatever the 
precise definition of rj, it is clear that the time interval between successive emitted 
Hawking quanta, is on average, large compared to the natural timescale associated 
with the energy of the individual emitted quanta. 

We now compare this with numerical estimates along the lines of Page’s results from 
the mid 1970’s [1-4]. The specific numbers will change, but the qualitative behaviour 
stays the same. Slightly modifying Page’s 1976 analysis [1, 2], for emission of a spin s 
quantum we define 

r = I T sim (u) (n) u dca; (4.7) 

(u>) = J uj T s£m (uj ) (n) u dw/T; (4.8) 

and 

—-— = 2 7r V] f T stm{u) {n) u —. (4.9) 

^7 binned J ^ 

Where available, we have used and adapted Page’s 1976 numerical results. Where not 
otherwise available, we have numerically estimated the greybody factors using product 
integral techniques adapted to the Regge-Wheeler/Zerilli potentials [26], followed by a 
numerical integration over frequencies. Our results are summarized in Table 1. 

Table 1. Semi-analytic and numerical estimates of g for massless bosons emitted from a 
Schwarzschild black hole. The semi-analytic estimates are from the current discussion. The 
numerical estimates are from Page’s 1976 results and our own extrapolations thereof. 


BOSONIC 

g 

V peak number 

V peak energy 

^/average frequency 

V binned 

Semi-analytic 

1 

327r 2 (2+W(—2e“ 2 )) 

327r 2 (3+W(—3e" 3 )) 

167T 6 

128 

27C(3) 

27f(3) 

405C(3) 2 

9 

Value 

1 

15.508 

27.465 

26.285 

14.222 

s = 0 

1 

20.65 

27.83 

26.78 

16.31 

s = 1 

2 

246.1 

259.1 

244.5 

216.3 

s = 2 

2 

5076 

5219 

4964 

4692 


Overall we see that introducing the greybody factors, (because they suppress T), 
always have the effect of increasing g. In particular, as one goes to higher spin there is a 
larger “angular momentum barrier” to overcome and the numerically estimated values 
of r] greatly exceed those obtained from the semi-analytic estimates based on a purely 
Planckian spectrum. For spin zero the semi-analytic estimates are in good agreement 
with the numerical results. For higher spin the semi-analytic estimates provide a lower 
bound, but generally sparsity of the Hawking flux is even more extreme than suggested 
by the semi-analytic estimate. 



































4.1.2 Massless Fermions 


The modifications for massless fermions are straightforward. First we note that the 
differential emission rate changes to 


dr 


g ck 2 

167T 2 exp (hck/ks) + 1 


dA; dA 


(4.10) 


(Note that the original version of the standard model of particle physics contained only 
massless chiral neutrinos with g — 1; as soon as one extends the standard model to 
contain non-chiral neutrinos then g — 2. We shall explicitly retain g, both for historical 
reasons and to facilitate comparisons with bosonic and Boltzmann results.) Integrating 
over wavenumber and area 

3gC(3) kjT 3 

167T 2 W’C 2 

For fixed g this is certainly less than the bosonic result, by a factor 3/4. Furthermore 
the peak in the number spectrum shifts upwards and is now at 



^ peak number 


ck 


peak number 


k B T 

h 


(2 + W (2e -2 )) , 


(4.12) 


while 


^ peak energy 


ck 


peak energy 


k B T 

h 


(3 + W (3e -3 )) . 


(4.13) 


Including both effects, for a Schwarzschild black hole, (where, as for the bosonic case, 
A = 27-nr 2 H and k B T H = hc/(4:7rr H )), we have 


V peak number 


128vr 2 (2 + W(+2e~ 2 )) 

81<7C(3) 


28.77434355... 

-> 1, 

9 


(4.14) 


and 


128vr 2 (3 + W(+3e~ 3 )) _ 40.62426089... ^ 

9 peak enero v XD 777X 1 • 


cenergy 81^(3) 

Using the average frequency 


9 


(uj) ->■ 


7vr 4 k B T H 


and so 


hi 


180 C(3) h ’ 

224tt 6 40.88835782... 


9 


I average energy 3645/'(3) 2 

Using the binned g yields 

256 28.44444444.. 


» 1. 


V binned 


9 g 


9 


> 1. 


(4.15) 


(4.16) 


(4.17) 


(4.18) 
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Thus it is clear that the time interval between emitted (massless fermionic) Hawking 
quanta is, on average, large compared to the natural timescale associated with the 
energy of the emitted quanta. We again compare this, where we can, with numerical 
results including the effects of greybody factors. See Table 2. Overall we see that the 
semi-analytic estimate is reasonably good for spin 1 / 2 , though we again expect that 
for higher spin the angular momentum barrier will depress T and enhance 77 . 

Table 2. Semi-analytic and numerical estimates of g for massless fermions emitted from a 
Schwarzschild black hole. The semi-analytic estimates are from the current discussion. The 
numerical estimates are from Page’s 1976 results and our own extrapolations thereof. 


FERMIONIC 

g 

V peak number 

V peak energy 

^average frequency 

V binned 

Semi-analytic 

g 

1287T 2 (2+W(2e -2 )) 

1287T 2 (3+1V (3e -3 )) 

224tt 6 

256 

81gC(3) 

81gC(3) 

3645gC(3) 2 

9.9 

Value 

g 

28.773 /g 

40.624/g 

40.88/(7 

28.444/(7 

Value 

2 

14.387 

20.312 

20.444 

14.222 

8 = 1/2 

2 

— 

29.3 

27.6 

— 


4.1.3 Massless Boltzmann particles 


For completeness, and for subsequent use, we now consider the case of particles which 
satisfy Boltzmann statistics. The differential emission rate changes to 

dT = -9— ck 2 exp ( — hck/k B T) d k dA. (4.19) 

8ir 2 


A key observation is 

1 ^ i 1 

gX _ 1 ^ e x ' > gX 1 

This implies (for fixed g ) that T bosonic > T bo itzmann > T fermionic- 
over wavenumber and area 


(4.20) 


In fact, integrating 


9 


k 3 B T 3 


T = ^~ -f— A, 


(4.21) 


Air 2 h 3 c 2 

which indeed lies between the bosonic and fermionic results. Furthermore the peak in 
the number spectrum also lies between the bosonic and fermionic results and is now at 


and 


peak number ck peak number 


2k B T t 


H 


h ’ 


peak energy ck peak energy 


3 k B T B 


h 


(4.22) 

(4.23) 
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In this case the ordering is cn peakbosonic < a7 peakbo itzmann < w peak fermionic implying 


bosonic ^ ?/boltzmann ^ 9 form ionic ■ 


(4.24) 


After a bit of work, for a Schwarzschild black hole we have 

_ 64vr 2 23.39461784. 

?? peak number 


V peak energy 


27 g 
32vr 2 


9 

35.09192677.. 


9# 9 

For the Boltzmann spectrum, the average energy is at the peak energy, so 

32t r 2 35.09192677... 


V 


average energy 


9 9 


9 


(4.25) 

(4.26) 


(4.27) 


Finally, the binned value of ?y in this situation agrees with that obtained from the peak 
in the number spectrum 


V binned 


647T 2 


23.39461784... 

9 


(4.28) 


We summarize our results in Table 3. Overall we see that the semi-analytic estimates 
for the sparsity of the Boltzmann flux fits nicely between those for bosonic and fermionic 
fluxes. 


Table 3. Semi-analytic estimates of ry for massless Boltzmann particles emitted from a 
Schwarzschild black hole. These semi-analytic estimates are from the current discussion. 


BOLTZMANN 

g 

^7 peak number 

V peak energy 

^average frequency 

V binned 

Semi-analytic 

g 

64:7V Z 

27 q 

32tt^ 

9.<? 

32tt^ 

9.9 

647 T z 

27 q 

Value 

g 

23.395 /g 

35.092 /g 

35.092/3 

23.395 /g 


4.1.4 The situation so far 

Up to this point we have only considered the emission of massless Hawking quanta from 
Schwarzschild black holes. We have seen that the semi-analytic estimate (based on 
assuming an exactly Planckian spectrum for the emission) gives a reasonably accurate 
slight under-estimate for the sparsity of s = 0 and s = I emission, but that for higher 
spin the angular momentum barrier quickly drives the sparsity parameter r] even higher. 
For photons numerical estimates including the effect of greybody factors yield rj > 200, 
for gravitons rj > 4500. 
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4.2 Reissner—Nordstrom black holes (uncharged quanta) 


In all the situations we have considered so far we have r/ oc 1/(T^Ah)', and we have 
then used properties of the Schwarzschild black hole to evaluate this quantity. For the 
more general Reissner-Nordstrom black holes, as long as we work in terms of the radius 
of the inner and outer horizons, r±, we can generically write 


r + — r_ 
= 2r\ ’ 


k bTh = 


Hck 

~2n' 


Ah = Anr 2 ^. 


(4.29) 


Consequently for the total emission of massless quanta, (which are then automatically 
guaranteed to be electrically neutral), and when working within the blackbody approx¬ 
imation, we have: 




.eissner- 


-Nordstrom ^Schwarzschild ^ 


(r+ - r_y 


T ^Schwarzschild* 


(4.30) 


The key point here is that adding charge to the black hole serves only to make the 
(semi-analytic estimate for the) Hawking flux even more sparse. 


4.3 Dirty black holes 


So-called “dirty black holes” are black holes surrounded by some generic matter fields [27]. 
Let us keep everything static and spherically symmetric and write 

dr 2 


ds 2 = —e 2< ^ r )(l — 2m(r)/r)c 2 d£ 2 + 


+ r\d0 2 + sin 2 # d<^ 2 ). (4.31) 


1 — 2 m(r)/r 

It is then straightforward to calculate the surface gravity and determine [27] 


k = 


e -<t>H 

2 r H 


( 1 - 8 ^rl) : Ah = 4 , Tr 


H- 


(Here pu is the energy density at the horizon.) Thence 


Ah 


V V Schwarzschild X 

But via the Einstein equations [27] 

47 tg n r 


(pH = 


- 87rG N p H r 2 H /c 4 ' 
(p - p r ) r 


'r H 


1 — 2 m(r)/r 


dr. 


(4.32) 


(4.33) 


(4.34) 


Now as long as the null energy condition [NEC] holds in the radial direction, then this 
quantity is guaranteed non-negative. Furthermore as long as the weak energy condition 
[WEC] is satisfied at the horizon, then we have pu > 0. Subject to these two classical 
energy conditions holding [28, 29] we have r) dirty black hole > ?? Schwarzschild- The key point 
here is that adding extra matter to the black hole environment serves only to make the 
(semi-analytic estimate for the) Hawking flux even more sparse. 
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4.4 Particle rest mass 


Adding particle rest masses for the emitted quanta one now has 

ck 2 


dr = 


9 


871-2 exp + c 2 k 2 jk B Tu^ =F 1 


d k dA, 


(4.35) 


where ujq = m^c 2 jfi is the Compton frequency. Unfortunately the relevant integral for 
the total number flux is no longer explicitly evaluable, at least not in any convenient 
form, and the best one can easily say is this: 

a k 3 T 3 

r = A h x |Li 3 (±l)| x f(hw 0 /k B T), (4.36) 

where Li n (z) denotes the polylog function, and 5 


/(*) 


roc x 2 

•'O exp(\/2 2 +x 2 )=Fl 


dx 


fc 


OC x 2 

0 exp(:r)=Fl 


dx 


< 1 . 


(4.37) 


Thus T certainly decreases. 

Similarly the location of the flux peak is no longer explicitly evaluable in closed 
form, (not even with the aid of the Lambert 1U-function), although it is possible to 
deduce on quite general grounds that the location of the flux peak moves upwards. 
These two observations are enough, however, to guarantee that 77 increases as one 
introduces particle rest masses, and so the semi-analytic estimate for the Hawking flux 
again becomes sparser. 

In counterpoint, the binned version of rj becomes 


1 


9 binned 


27T 


i dr 

a/wq + c 2 k 2 dk 


dk, 


(4.38) 


whence, more explicitly 6 


1 _ g_ [ 1 _ ck 2 A h _ 

9 binned 47T J y/ Uq + C 2 k 2 eX p + C 2 k 2 /k B T H ^ =F 1 


(4.39) 


It is manifest from the above that I /77 binned again decreases, so we see that (the semi- 
analytic estimate for) //binned again increases, as particle rest masses are added to the 
emitted Hawking quanta. However, as we shall now see, once one enters the super¬ 
radiant regime the discussion becomes much messier. 

5 An explicit representation of the integral that appears in the numerator of f(z) can be obtained 
in terms of an infinite sum of Bessel functions, but that expression does not seem to be particularly 
useful. 

6 An explicit representation of this integral is also obtainable in terms of an infinite sum of Bessel 
functions, but this expression does not seem to be particularly useful. 
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4.5 Summary 


At least for spherical symmetry, and in the non-super-radiant modes, the semi-analytic 
estimates for the sparsity of Hawking emission seem to be bounded below by that for 
the Schwarzschild black hole. Insofar as we have been able to check the numerical 
estimates based on including the effects of greybody factors, this continues to hold in 
more general situations; the Schwarzschild black hole provides a good lower bound on 
sparsity, at least in the absence of super-radiance. 


5 Super-radiant situations 


Super radiance can occur when for one reason or another one has a chemical potential 
/i. The bosonic occupation number then becomes 


(n) u, -)• 


1 

exp {(ftw - n)/k B T H } 


1 


(5.1) 


If the chemical potential becomes positive (/i > 0) then the occupation number diverges 
for ftw = /i and becomes negative for to e (0, p/h). In all situations of interest the 
greybody factor Te sm (co) simultaneously exhibits a zero and then becomes negative 
so that the product Tz srn (ijj) (n) ^ remains finite and positive [1]. The region where 
both occupation number and greybody factor are negative is called the “super-radiant 
regime”. This phenomenon occurs only for bosons, and (in the current black hole 
context) arises either due to charged particle emission from a charged black hole, or 
for low L z angular momentum modes from a rotating black hole (or a combination of 
these two effects). For a more general discussion, not necessarily black-hole related, see 
for instance [30], and references therein. 

There is a definitional ambiguity as to whether one should include the super-radiant 
modes as part of the Hawking flux, or treat them separately. For instance, super¬ 
radiance has a well-known classical limit — the Penrose process for mining black hole 
rotational kinetic energy. If one insists on considering the super-radiant modes as a 
subset of Hawking radiation, then one would also be forced to assert that the Penrose 
process is the classical limit of the Hawking process — an assertion that many would 
feel is excessive in its generality. 

Certainly both the super-radiant modes and non-super-radiant modes are each 
subsets of “quantum vacuum radiation”, but the underlying physics is very different 
and we feel there are good reasons to keep the two notions distinct. Hawking quanta 
are related to horizons, while super-radiance is related to (generalized) ergo-regions. 
Hawking quanta are temperature dominated, while super-radiant quanta are chemical 
potential dominated. We shall reserve the phrase “Hawking modes” for the non-super- 
radiant subset of the quantum vacuum radiation modes. 
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5.1 Reissner—Nordstrom black holes (charged quanta) 

The key change for the emission of charged quanta from a charged black hole is that the 
boseonic occupation factor picks up a “chemical potential” /i = qV E depending on the 
charge q of the emitted quanta and the electromagnetic potential Vh at the horizon: 


exp {(hw - qV H )/k B T H } - 1 


(5.2) 


If the mass of the emitted particle satisfies m^c 2 < qV B then there is a critical frequency 
at which the bosonic occupation number diverges (and below which the bosonic occu¬ 
pation number is actually negative). In that super-radiant region there is no escaping 
the relevance of the grey-body factors — the T( mn (u) non-perturbatively differ from 
unity, and there is no longer any good reason to really trust the semi-analytic approx¬ 
imation based on a pure Planck spectrum. In this context, we can in fact extend the 
usual notion of a ergo-surface in a particle-dependent manner by defining the radius 
t'e of the ergo-surface to be 

m 0 c 2 = qV(r E ), (5.3) 

while the associated ergo-region is defined by 

{r : m 0 c 2 < qV(r)}. (5.4) 


Note that these generalized ergo-regions depend on both rest mass and electric charge. 
Astrophysically these generalized ergo-regions are of limited usefulness, (since astro- 
physical black holes tend to quickly neutralise). Accordingly, we find it most useful to 
concentrate the discussion on the Kerr spacetime. 


5.2 Kerr and Kerr—Newman black holes 


For the Schwarzschild black hole in the pure Planck blackbody approximation we saw 
that rj oc 1 /(T^jA h )\ and we then used properties of the Schwarzschild black hole to 
evaluate this quantity. For the more general Kerr and Kerr-Newman black holes, as 
long as we work in terms of the radius of the inner and outer horizons, r±, we can 
generically write 


r + — r_ 
2 (r 2 + a 2 ) ’ 


k B Tn — 


hen 


A h = 4n(r 2 + + a 2 ). 


(5.5) 


Therefore, within the context of the semi-analytic approximation 


(ft Ah) Kerr — Ah) Schwarzschild X 


(r+ — r_) 2 
r 2 + a 2 


(5.6) 
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Consequently for the total emission of massless quanta, (which are then automatically 
guaranteed to be electrically neutral), and when working within the blackbody approx¬ 
imation, we have: 


h Kerr V Schwarzschild X . .. T Schwarzschild • 

1 T , - T I ^ 



(5.7) 


(r+ — r_) 


The key point here is that adding charge and/or angular momentum to the black hole 
serves only to make the (semi-analytic estimate for the) Hawking flux even more sparse. 

In contrast, Hod [14] asserts the equivalent of r\ = 0(1) for highly extremal Kerr. 
(That is, for a = J/M 2 = a/M < 1.) Hod achieves this by lumping super-radiance 
into the Hawking flux. In the extremal limit (k —>■ 0) the greybody factors are approx¬ 
imately [ 1 ] 


T(uj) & C\ s (A h uj[uj - 


(5.8) 


The bosonic occupation number is 


(5.9) 


exp {h(uj — mfln) I^bTh} — 1 


These both change sign at uj = 

• For 0 < uj < the greybody factor is negative. 

• For 0 < uj < rniljr the bosonic occupation number is negative. 

Consequently, for 0 < oj < mO,H the super-radiant emission is not well-approximated 
by a blackbody. So it makes sense to split the integral into two regions: 

• 0 < oj < mil h — super-radiant emission. 

• mil h < to < oo — Hawking emission. 

It is well-known that in the extremal limit (k —)■ 0) super-radiance is known to dominate 
over the Hawking flux [1-5]. If we concentrate on the binned version of r/ (which 
for present purposes is more important than the flux) then slightly modifying Page’s 
analysis 



£m 


(5.10) 


Then in the near-extremal limit: 



(5.11) 
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Here: x = w/(mfijj) and e = {hmVL h) / (k bT h ) 1 . As usual, the emission is 

dominated by the lowest available angular momentum state £ = m = s: 


- « (A h s 2 Q 2 h ) 2s+1 
V 


(x[x — 1]) 


2s+l 


dx 


J exp(e[x — 1]) — 1 x 
Splitting into super-radiant and Hawking modes: ' 


and 


V super-radiant 

1 


c a ,s (A H s 2 n 2 H ) 


2q2 \2s+1 


V Hawking 


c s ,s (A H s 2 n 2 H ) 


2q 2 \2s+l 


1,1 (x[x — l]) 2s+1 dx 
Jo exp(e[x - 1]) - 1 x ’ 

f°° (x[x — l]) 2s+1 dx 
/, exp(e[a: — 1]) — 1 x 


(5.12) 


Recalling that e 1, some brute-force integration then yields the estimates: 

super-radiant = O (c°) ; I(H>wking = O (f 2 ‘ +2 ) »» 1. 


(5.13) 


(5.14) 


(5.15) 


So super-radiance dominates in the extremal limit. In fact, super-radiance leads to 
rapid spin-down with small energy loss [1-5], until the system goes non-super-radiant, 
and then the “normal” Hawking effect takes over. The quantitative details are messy, 
but the overall message is clear: Sparsity of the Hawking flux is the dominant feature 
of the Hawking evaporation process. 


6 Discussion 

So long as the temperature remains below the rest energy of the electron we can safely 
ignore emission into charged particles. Indeed, considering emission into individual 
angular momentum modes makes sense only within the specific framework of the Page 
analysis [1-3], and such complications can be safely ignored as long as the Hawking 
radiation can be treated as an approximate blackbody. Collecting all the various results 
we see that throughout the initial low-temperature stages of Hawking evaporation: 

• For bosonic Hawking quanta rj 1 in all situations of interest. 

• For fermionic Hawking quanta rj 3> 1 in all situations of interest. 

7 Note the split takes place “in parallel”; that is 

1 1 1 

— -- 1 -. 

V V super-radiant Hawking 
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• For Boltzmann Hawking quanta r/ 1 in all situations of interest. 


This is best summarized as the statement that the Hawking flux is extremely sparse — 
the average time between emission of Hawking quanta is very large compared to the 
timescale set by the energies of the Hawking quanta. The Hawking quanta are dribbling 
out one-by-one, with very large interstitial gaps. This phenomenon persists throughout 
the entire evolution of the black hole, both in early stages and (with suitable minor 
modifications) in late stages. (Earlier authors have largely confined their attention to 
the late-time high-temperature regime [ 6 - 11 ].) 

The bad news is this: Compared to numerics the semi-analytic estimates are often 
off (under-estimating sparsity) by factors of 100 or even more. Sparsity if anything 
increases when you do numerics that includes the effects of greybody factors. Why 
the discrepancy? Individual photons can at best be localized (both in the direction 
of motion and transversely) to about a wavelength. The usual blackbody emission 
estimate is based on treating the photons as particles, not waves. Once the physical 
size of the emitter becomes smaller than a “typical wavelength”, (more exactly, the 
thermal wavelength), the usual blackbody emission estimate is no longer trustworthy, 
providing at best a qualitative estimate, typically a lower bound on sparsity. But note 
that the usual statistical mechanics result 


dT 


g ck 2 

87 r 2 exp( hck/ksTn ) — 1 


d k dA. 


( 6 . 1 ) 


should still be valid for u —>■ 00 . So semi-analytic estimates are OK in the Boltzmann 
tail, but may mis-estimate things in the low-frequency regime. Sparsity, however, is 
here to stay — modulo arguments on how to classify super-radiance. 


The sparseness of the Hawking flux has a number of perhaps unexpected kinematical 
implications: 

• While early-stage Hawking radiation from Schwarzschild or Reissner-Nordstrom 
black holes is spherically symmetric, this spherical symmetry is only a long-term 
statistical statement obtained after averaging over very many Hawking quanta. 

• Early-stage Hawking evaporation should be seen as a long chain of independent 
2-body decay processes involving photons, gravitons, and neutrinos. (Similarly, 
late-stage Hawking evaporation, once the temperature exceeds Aqcd, should be 
viewed as a long chain of 2-body decay processes proceeding via the emission of 
hadronic jets [6-11].) 
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• When analyzing the emission of individual Hawking quanta one should use the 
special relativistic kinematics that is applicable in the asymptotic spatial region. 
Consider a black hole initially of rest energy rriiC 2 which emits (in its rest frame, 
but as seen from spatial infinity) a particle of energy hu and rest energy Huq, 
thereby reducing its rest energy to rrifC 2 . Then the conservation of 4-momentum, 
(when written in terms of the 4-velocities V), implies 

(me 2 ) Vi = {m f c 2 ) V f + (hu 0 ) V H . (6.2) 


Therefore 

(mjC 2 ) 2 + (hu q) 2 — 2 (miC 2 )(hu) = (m/C 2 ) 2 , (6.3) 


and so 


, 2\2 


hu = 


(rriiC 2 )~ + ( hu 0 ) 2 - ( m f c 2 ) 
2 rrijC 2 


(6.4) 


Depending on whether or not one views black hole masses as being quantized or 
continuous, one can view this either as a normal 2-body decay, or as the decay 
of one IMP (“indefinite mass particle” ~ unparticle) into another IMP [31, 32]. 
It may be profitable to reconsider and reanalyze the entire Hawking evaporation 
process from this point of view. 8 


• More radically, the fact that the Hawking flux is so sparse calls into question 
the appropriateness of the use of bosonic and fermionic statistics, (at least for 
non-equilibrium Hawking evaporation into the Unruh vacuum state). For the 
intrinsically equilibrium Hartle-Hawking vacuum there is absolutely no doubt 
that bosonic and fermionic statistics are the relevant ones, but when quanta 
are so well separated as they are in the non-equilibrium Hawking evaporation 
process it is less clear that bosonic and fermionic statistics make sense, since 
these normally require multiple occupation of the same quantum mode, typically 
achieved by being at (or at least temporarily near) equilibrium. 9 On the other 
hand, in the super-radiant regime it is clear that bosonic statistics is essential , 
and that the blackbody approximation (Planck approximation) for the emitted 
flux is a poor one. 

8 MV’s early views on the importance of the Hawking cascade can be found in reference [33]. 
Although MV is no longer in favour of the particular way that black hole entropy was discretized in 
that article, the comments regarding the importance of the final “particle cascade” leading to complete 
evaporation of Planck-scale black holes still hold. 

9 For instance, this certainly happens in a “quasi-continuous-wave” laser, where photon emission 
is not really continuous — finite bandwidth is inversely related to finite coherence time. Typically 
one has short microsecond photon pulses with very high occupation numbers, so bosonic statistics is 
completely appropriate. 
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• Finally, and even more speculatively, the sparsity of the Hawking flux provides a 
very obvious place to hide information without significant energy cost — the gaps 
between the Hawking quanta can easily encode significant information without 
disturbing the (time-averaged) emission spectrum — this would effectively be a 
“phase modulated” Hawking flux. 

In summary: The extreme sparsity of the Hawking flux is a rather under-appreciated 
feature of the Hawking evaporation process. Hawking evaporation is seen to be a slow 
dribbling out of Hawking quanta, with the individual quanta being well separated in 
time. (The gap between individual quanta being hundreds of times longer than the 
natural timescales associated with the quanta.) This “decay chain” viewpoint has a 
number of tantalisingly interesting physical implications. 
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